A Pod Projection Method for Large-scale Algebraic Riccati Equations

The solution of large-scale matrix algebraic Riccati equations is important for instance in control design and model reduction and remains an active area of research. We consider a class of matrix algebraic Riccati equations (AREs) arising from a linear system along with a weighted inner product. This problem class often arises from a spatial discretization of a partial differential equation system. We propose a projection method to obtain low rank solutions of AREs based on simulations of linear systems coupled with proper orthogonal decomposition. The method can take advantage of existing (black box) simulation code to generate the projection matrices. We also develop new weighted norm residual computations and error bounds. We present numerical results demonstrating that the proposed approach can produce highly accurate approximate solutions. We also briefly discuss making the proposed approach completely data-based so that one can use existing simulation codes without accessing system matrices.